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Abstract 

The gauge covariant lattice Dirac operator D which has recently been proposed by 
Neuberger satisfies the Ginsparg-Wilson relation and thus preserves chiral symmetry. 
The operator also avoids a doubling of fermion species, but its locality properties 
are not obvious. We now prove that D is local (with exponentially decaying tails) 
if the gauge field is sufficiently smooth at the scale of the cutoff. Further analytic 
and numerical studies moreover suggest that the locality of the operator is in fact 
guaranteed under far more general conditions. 
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1. Introduction 


Many technical complications in the standard formulation of lattice QCD have to do 
with the fact that chiral symmetry is violated at the scale of the cutoff. In particular, 
the quark masses are not protected from additive renormalizations and the leading 
lattice effects in physical amplitudes are proportional to the lattice spacing a rather 
than being of order 

Somewhat surprisingly it has recently turned out [2-6] that chiral symmetry can 
be preserved on the lattice, without fermion doubling, if the lattice Dirac operator 
D satisfies a certain algebraic relation, 

= aDj^D, ( 1 . 1 ) 

originally due to Ginsparg and Wilson [1]. We shall not discuss the significance and 
consequences of this identity here, but refer the reader to the original papers quoted 
above and the rapidly growing literature on the subject [7-18]. A point which should 
be emphasized however is that the Ginsparg-Wilson relation only guarantees that 
the lattice theory has the same chiral symmetries as the continuum theory. Locality, 
the correct behaviour in the classical continuum limit and the absence of doubler 
modes are additional constraints which any decent lattice Dirac operator should 
satisfy. 

Starting from the overlap formulation of chiral gauge theories, a relatively simple 
solution of the Ginsparg-Wilson relation has been found by Neuberger some time 
ago [4]. Explicitly it is given by 

D = -{l-A{A^A)-^/n, A = l + s-aD^, (1.2) 

a 

where denotes the standard Wilson-Dirac operator, 

Dw = + - aV*^V^} , (1.3) 

and s is a real parameter in the range jsj < 1 which will be fixed later (cf. appendix A 
for unexplained notations). Neuberger’s operator is manifestly gauge covariant and 
can be shown to have no doubler modes. One can also easily verify that it converges 
to the expected expression in the classical continuum limit, up to a finite normaliza¬ 
tion constant, but the requirement of locality is not obviously fulfilled. Evidently it 
is very important to check that the operator is local because the universality of the 
continuum limit depends on this fundamental property. 
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Before we begin with the detailed discussion of Neuberger’s operator it may be 
helpful to state what precisely is meant if we say that D is local. Strict locality 
would imply that the non-zero contributions to the sum 

D^{x) = D{x,y)i;{y) (1.4) 

y 

come from the points y in a finite neighbourhood of x. Moreover the kernel D{x, y) 
should only depend on the gauge field variables residing near x. From eqs. (1.2) and 
(1.3) it is obvious, however, that Neuberger’s operator is not local in this restricted 
sense. A more general definition of locality is hence adopted here, where the kernel 
is allowed to have exponentially decaying tails at large distances. As long as the 
rate of decay can be shown to be proportional to the cutoff 1/a, the sum in eq. (1.4) 
will be completely dominated by the contributions from a bounded region around x 
with a fixed diameter in lattice units. In particular, from the point of view of the 
continuum limit there is little doubt that this kind of locality is as good as strict 
locality. 

It is clear that a non-locality of the Neuberger operator can only arise from the 
inverse square root of A in eq. (1.2). Most of the time we shall thus be concerned 
with the properties of this operator. In section 2 we first consider the case where 
AtA is bounded from below by a positive constant. Expanding the inverse square 
root of A^^A in a series of Legendre polynomials, the locality of D may then be 
proved straightforwardly. Moreover, the required lower bound on the spectrum of 
AtA can be established rigorously if the gauge field is sufficiently smooth at the 
scale of the cutoff and Neuberger’s operator is hence guaranteed to be local for all 
these fields. 

One might expect that D becomes increasingly non-local when A^A develops a 
zero mode. We briefly examine this question in section 2, using series expansions, 
and find that this is actually not so in general. The numerical studies reported in 
section 3 confirm this and they also provide a realistic estimate of the localization 
range of the operator at the gauge couplings of interest. All these results fit into a 
simple picture which suggests that Neuberger’s operator (with an appropriate choice 
of the parameter s) is local for all statistically relevant gauge field configurations. 
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2. Spectrum of A and locality of D 


In this section some rigorous results are established which show that D is local (in 
the sense explained above) under certain conditions. Along the way we shall find 
that the locality of D is closely related to the spectral properties of A. This leads 
to interesting qualitative insights which allow us to be more confident about the 
interpretation of our numerical studies. 

2.1 Series expansion o/(A^A)“^/^ 

Following a suggestion of Bunk [19] the inverse square root of A^ A may be expanded 
in a series of Legendre polynomials. The expansion is suitable for numerical appli¬ 
cation, but here we use it as a theoretical tool to discuss the locality properties of 
the Neuberger operator. In the following lines the detailed form of A^A does not 
matter. To ensure the convergence of the Legendre expansion we however assume 
that the bounds f 

u<A’fA<u (2.1) 

hold for some strictly positive constants u < v. Whether this is the case for a given 
gauge field is a separate issue which will be addressed later. 

The Legendre polynomials Pk (z) may be defined through the generating function 

OO 

+ = J2t’"Pk{z). ( 2 . 2 ) 

fc =0 

Usually z is taken to be a number, but eq. (2.2) remains meaningful if we substitute 

z = {v + u — 2A'^ A)/{v — u). (2.3) 

It is easy to show that this operator has norm less than or equal to 1. In par¬ 
ticular, using the operator calculus and the well-known properties of the Legendre 
polynomials (as quoted in ref. [20], for example), this implies 

ll^fc(^)ll<l- (2.4) 

The expansion (2.2) is hence norm convergent for all t satisfying |t| < 1. 

I Here and below an inequality between operators stands for the corresponding inequality between 
the expectation values of the operators in arbitrary normalizable states 
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We now introduce a parameter 9 through 


cosh.9 = {v + u)/{v — u), 6* > 0, (2.5) 

and set t = e“®. Eq. (2.2) then assumes the form 

OO 

= K Y.t^Pkiz), K = {4t/{v - u)}'/" , (2.6) 

k=0 

since 1 — 2tz +is proportional to A'^A for this choice of t. 

2.2 Legendre expansion and locality of D 

We now show that the convergence of the Legendre expansion (2.6) and the strict 
locality of A imply that D is local (with exponentially decaying tails). The lattice is 
here taken to be infinitely extended in all directions. In view of eq. (1.2) it suffices 
to establish the locality of the inverse square root of A"^A. 

The kernel G{x, y) which is associated with this operator, 

{A^ A)-^/‘^il){x) = a^Y^ G{x, y)i>iy), (2.7) 

y 

is a matrix acting on the Dirac and colour indices of the fermion field at the point y. 
If we define the kernels Gkix,y) representing the operators Pkiz) in the same way, 
we have 

OO 

G{x, y) = t’^Gkix, y). (2.8) 

k=0 

It is easy to show that the norm convergence of the Legendre expansion implies the 
absolute convergence of this series for all points x and y. Actually, from eq. (2.4) 
one infers that 

||Gfc(x,y)|| < 1 for all k,x,y, (2.9) 

where the norm here is the matrix norm in Dirac and colour space. 

We now note that Gk{x,y) vanishes unless x and y are sufficiently close to each 
other, because Pk{z) is a polynomial in A^A and A is a combination of nearest- 
neighbour difference operators. If we introduce the “taxi driver distance” 

\\x-y\\i = '^\x^j,-y^,\, ( 2 . 10 ) 

n 


4 



the precise statement is that 


Gk{x,y)=0 for all /c < ||x — /2a. (2-11) 

Restricting the sum (2.8) to the non-zero terms and recalling eq. (2.9), the bound 

a^\\G{x,y)\\ < Y^exp{-6l||x-?/||^/2a} (2.12) 

is thus obtained. In particular, the kernel is exponentially decaying at large distances 
with a rate proportional to the cutoff 1/a. Moreover its dependence on the gauge 
field is local in a similar way, i.e. up to exponentially small tails. 

A technical detail we wish to emphasize is that the localization range 2a/9 and 
the proportionality constant in eq. (2.12) only depend on the bounds u and v. As 
long as these can be chosen uniformly in the gauge field, the Neuberger operator is 
guaranteed to behave essentially as a strictly local operator. The differentiability of 
the operator with respect to the gauge field can also be proved under these conditions 
(appendix B). 

2.3 Bounds on A^A 

The proof of the locality of D given above depends on the convergence of the Legen¬ 
dre expansion and thus on the validity of the bounds (2.1) for some positive constants 
u < V. The upper bound is easily seen to hold for any gauge field if we rewrite A in 
the form 

A = 1 -I- s -I- ^ {i(l - 7 ^)aV^ - |(1 -I- 7 ^)aV^} . (2-13) 

Using the triangle inequality it then follows that 

||A|| = ||.4t||<8 (2.14) 

for |s| < 1 and A^A is hence uniformly bounded from above. 

As far as the lower bound is concerned it is clear that A^ A is non-negative, but one 
knows that the operator can have zero modes for some gauge field configurations. A 
uniform lower bound is hence excluded. The strict positivity of A^ A may however 
be established if the gauge field is sufficiently smooth at the scale of the cutoff. To 
make this more precise, let us suppose that 

||1 — C/(p)|| < e for all plaquettes p, (2-15) 
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where U (p) denotes the product of the gauge field variables around p and the norm 
is the matrix norm in colour space. For s = 0 the inequality 

y4^y4>l-30e (2.16) 

may then be proved (appendix C) and a similar bound may be deduced from this 
for general s by substituting ^4 = 5 + A\s-q and using triangle inequalities. A is 
hence uniformly bounded from below by a positive constant if e is sufficiently small. 
In particular, the locality of D is guaranteed under these conditions. 

2-4 Locality of D in the presence of near-zero modes of A^ A 

When the gauge field strength is not small it can happen that some of the eigenvalues 
of A^A are very close to zero or even equal to zero. Since the exponent 6 defined in 
subsect. 2.1 is proportional to at small u, one is tempted to conclude that the 
locality of D is lost in this situation. We now show that this is not so in general. 

Let us consider the case where the spectrum of A^A is contained an interval [u, u] 
except for an isolated eigenvalue A in the range 

0 < A < |u. (2.17) 

As before the constants u < v are some fixed positive numbers while A is allowed to 
become arbitrarily small. The projector on the associated eigenspace is given by 



(2.18) 


where the integration contour is a circle in the complex plane centred at the origin 
with radius |tt. We now first prove that P is local by noting that the operator in 
the square bracket on the right-hand side of the identity 

(u; - A^A)-^ = {w* - A^A) [(«;* - A'^A){w - A^A)] (2.19) 

has eigenvalues between {\u)^ and {v juY. Its inverse may thus be expanded in 
a rapidly convergent series of Chebyshev polynomials and it follows from this that 
the kernel associated with P is local with exponentially decaying tails. It should be 
emphasized that this remains true even if A approaches zero, which is a regular case 
in the above equations. In particular, the localization range of P is determined by 
u and V alone. 
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To establish the locality of D we write 


(^t^)-i/2 ^ + (ylt^)-i/2(i _ p) (2.20) 

and expand the terms on the right-hand side of this equation in Legendre polyno¬ 
mials. In particular, 


{A^A)-^/\l -P)=K Y,t'^Pk{z){l - P), (2.21) 

fc =0 

where t, z and k are as given in subsect. 2.1. The convergence of the series is 
guaranteed since ||Pfc( 2 ;)(l — T’)|| < 1. Proceeding essentially as in subsect. 2.2, this 
implies the locality of the operator (2.21). For the Legendre expansion of the hrst 
term in eq. (2.20), the spectral bounds u,v should be replaced by u,v, where 

V = 2u, u < X < V. (2.22) 

The associated expansion parameter t is independent of u and the localization range 
of the operator is hence the same for all values of A. Note that the divergence of 
the normalization factor k cancels when the inverse square root of A^ A is multiplied 
with A. In particular, Neuberger’s operator remains hnite and local in the limit 
A ^ 0. 

2.5 Summary 

Perhaps the most important result obtained in this section concerns the small held 
region where the bound (2.15) holds for some e strictly less than Neuberger’s 
operator is local in this case with exponentially decaying tails. Moreover the lo¬ 
calization range is uniformly bounded from above by a constant depending on e 
only. 

In the large held region the situation appears to be more complicated and it could 
be that the locality of D cannot be guaranteed for all helds. Nevertheless we have 
been able to show that the presence of near-zero modes of A^A does not by itself 
imply any non-locality. Our analytical investigations rather suggest that D remains 
local as long as the continuous spectrum of A"^ A is separated from zero by a positive 
gap. 
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3. Numerical studies of Neuberger’s operator 


In numerical simulations of quenched lattice QCD, using the Wilson plaquette ac¬ 
tion, the representative gauge field configurations at the gauge couplings of interest 
have relatively large average plaquette values. As a consequence the locality of D 
cannot be guaranteed on the basis of the results of section 2 alone. The purpose of 
the numerical studies reported here is to obtain some direct evidence for (or against) 
the locality of D in this situation and to check whether the qualitative picture is as 
suggested by the theoretical analysis. 

The lattices that we have considered are of size L in all directions with periodic 
boundary conditions. L/a has been set to 12 or 16. The gauge group is taken to be 
SU(3), with three values of the bare coupling = 6//3 corresponding to /3 = 6.0, 6.2 
and 6.4. Following standard procedures, a representative ensemble of statistically 
independent gauge field configurations has been generated for each lattice. It should 
be emphasized that all results refer to quenched QCD. In the full theory the situation 
could be different, although there is currently no reason to expect this. 


3.1 Localization range of D 
Let us consider the source field 


J 1 \i X = y and a = 1, 
\ 0 otherwise. 


(3.1) 


where y is some particular point on the lattice and a runs over the colour and Dirac 
indices of the field. We are then interested in the decay properties of 


ip{x) = A{A^A) ^/^r/(x) 


(3.2) 


at large distances \\x — y\\i [cf. eq. (2.10)]. It is implicitly understood here that the 
coordinate differences x^ — y^ are taken modulo L so as to minimize the distance. 
In particular, the largest possible distance is 2L. 

In fig. 1 we plot the expectation value of the function 

/(r) = max{||V’(x)|| | jjx - yjji = r] (3.3) 

for various values of s. The norm ||^/)(x)|| in this definition is the usual vector norm. 
To compute ^l>{x) we have used a Chebyshev approximation for the inverse square 
root of A'^A, with coefficients adjusted so that a relative accuracy better than 10“® 
is achieved [21,22]. One needs to know the extremal eigenvalues of A^A for this, but 



f(r)) 

1 

10-1 

10-2 

10-3 

10-4 

10-5 

10-6 

10 -^ 

10-6 


Fig. 1. Expectation value of f{r) [eq. (3.3)] as a function of the distance r. Open 
symbols correspond to s = 0 while the filled symbols represent the data at s = 0.4 
(/3 < 6.2) and s = 0.2 (/3 = 6.4) respectively. The statistical errors are not visible on 
the scale of this plot. 

as discussed below they can be calculated reliably with a modest effort. A technical 
point we wish to emphasize is that the relatively high numerical precision quoted 
above is required to avoid systematic effects in the calculated values of /(r) at large 
distances. 

In all cases considered {f{r)) is rapidly decaying when the distance r increases. 
Finite-volume effects appear to be negligible here and significant differences between 
the curves at different /3 and s are only seen when r/a is larger than 10 or so. For 
r/a > 13 the data can be represented by a single exponential, 

(/(r))oce—(3.4) 

with exponents u as listed in table 1. One may be worried at this point that the 
fluctuations of /(r) are large, but our experience is that the mean deviations are at 
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Table 1. Values of exponent v [eq. (3.4)] 


/3 

Lja 

s 

V 

6.0 

12 

0.0 

0.28 


12 

0.4 

0.49 


12 

0.6 

0.45 

6.2 

12 

0.0 

0.35 


12 

0.4 

0.49 


12 

0.6 

0.42 

6.4 

16 

0.0 

0.40 


12 

0.0 

0.40 


12 

0.2 

0.53 


12 

0.4 

0.49 


most a fraction of the average value and would thus be barely visible in fig. 1. 

From the figure and the table it is evident that (/(r)) becomes nearly independent 
of /3 if s is chosen appropriately. The lowest curve in fig. 1 is also practically matched 
by the curve that one obtains in the free quark theory. In other words, as far as 
(/(r)) is concerned, the localization properties of D at /3 > 6.0 and with a good 
choice of s are roughly the same as in the free case. 

That some tuning of s is required to preserve the localization range of D does 
not come as a total surprise, because the spectrum of the Wilson-Dirac operator 
moves to the right in the complex plane when /3 is decreased. In particular, the 
critical bare mass rric is shifted to —0.68/a, —0.74/a and —0.82/a at /3 = 6.4, 6.2 
and 6.0 respectively [23]. Since 

A =-a{D^ + mo), rriQ = -{I + s) / a, (3.5) 

a positive value of s partly compensates for this and ensures that rrio keeps away 
from nic by an appreciable margin. 

3.2 Spectrum of A 

To make contact with the theoretical discussion of section 2 we now proceed to ex¬ 
amine the distribution of the low-lying eigenvalues of A'^A. For any given gauge field 
configuration these eigenvalues can be computed by minimizing the Ritz functional 
using a conjugate gradient algorithm. The method has previously been applied and 
is described in detail in refs. [24,25]. 
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Fig. 2. Monte Carlo time evolution of the four lowest eigenvalues of A^A. The 
smallest eigenvalue is represented by an open square and the higher ones by filled 
squares. All data are from the 12^ lattice with s = 0. 

As shown by fig. 2 the spectrum of A^A depends quite strongly on the gauge 
coupling. At /3 = 6.4 the lower end of the spectrum appears to be clearly separated 
from zero and A thus satisfies the bounds (2.1). This is not so at /? = 6.2 and 6.0, 
where one has a non-zero probability to find eigenvalues that are orders of magnitude 
below the rest of the spectrum. Moreover the band of the ordinary modes is wider 
at these couplings. 

When s is set to higher values (such as those quoted in fig. 1) the qualitative 
features of the spectrum do not change, but the level of the ordinary low-lying 
eigenvalues is raised by a factor 2 or so. The maximal eigenvalue, on the other hand, 
hardly changes and stays around 41 for all (3 and s that we have considered. So far 
we have only been able to analyse a limited number of gauge field configurations 
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on the larger lattice and thus cannot make a detailed statement about the volume 
dependence of the spectrum at this point. There is, however, a clear tendency that 
the probability for near-zero modes increases with the lattice size. 

3.3 Localization properties of the near-zero modes 

The minimization of the Ritz functional not only yields the low-lying eigenvalues of 
A but also the corresponding eigenfunctions. In particular, the localization prop¬ 
erties of the near-zero modes can be studied straightforwardly. A possible definition 
of the localization range of a given wave function is discussed in ref. [27]. Proceeding 
along these lines we have found that all near-zero modes are well localized with ex¬ 
ponentially decaying tails. Moreover we have observed that the localization ranges 
shrink significantly when s is increased. All this completely agrees with the results 
previously reported by Edwards et al. [26] (see also ref. [27]). 

For illustration let us consider an eigenvector <p{x), suitably normalized, with 
maximal magnitude ||(/)(x)|| at some point x = y. The function 

h{r) = max {||(()(x)|| | ||x — 2/||i = r} (3-6) 

then provides an upper bound on the wave function at distance r from the centre of 
its localization region. A typical result for h{r) is plotted in fig. 3. The associated 
eigenvalue is nearly two orders of magnitude below the band of the ordinary low-lying 
modes in this example. 

In fig. 3 the function /(r), calculated for the same gauge field configuration and 
with the same choice of y, is also shown. Comparing with (/(r)) one clearly sees 
that the presence of the near-zero mode does not affect the localization properties 
of D in any significant way. 

3-4 Synthesis 

Taken together the theoretical and numerical results reported in this paper suggest 
that the spectrum of A^A is clearly separated from zero at large f5, for all statistically 
relevant gauge field configurations. The locality of Neuberger’s operator is rigorously 
guaranteed in this situation. It is our experience, however, that the theoretical 
bounds tend to over-estimate the localization range by a large factor. At f} = 6.4 
and s = 0, for example, the exponent 612 = 0.03 which one obtains by inserting 
the numerically determined spectral bounds in eq. (2.5) is much smaller than the 
exponent v = 0.4 quoted in table 1. 

At the lower values of (3, near-zero modes develop and become increasingly fre¬ 
quent, but they are well localized (if s is chosen appropriately) and thus do not 
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Fig. 3. Plot of the magnitude h{r) of a typical near-zero mode at /3 = 6.2, s = 
0.4 and Lja = 12, together with the function /(r) and the corresponding ensemble 
average {f(r)). 

destroy the locality of Neuberger’s operator. On the lattices that we have studied 
the near-zero modes are isolated from the rest of the spectrum. The observed lo¬ 
cality properties are, therefore, completely in line with the theoretical discussion of 
section 2, although here again the analytical estimates of the localization radius are 
far too pessimistic. 

An important point to note in this context is that the probability to find local¬ 
ized near-zero modes grows proportionally to the volume, because widely separated 
regions on the lattice basically behave as separate systems. Eventually there will be 
configurations with several near-zero modes and in the infinite volume limit these 
modes form a dense so-called pure point spectrum (cf. ref. [29], chapter VII). It seems 
unlikely that this statistical phenomenon has any relevance for the locality proper¬ 
ties of Neuberger’s operator, but it is clearly desirable to check this by extending 
our studies to larger lattices. 
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4. Concluding remarks 


At this point there is little doubt that Neuberger’s operator is local at small gauge 
couplings gQ, although a rigorous proof of the locality is only available if the gauge 
field is assumed to satisfy the smoothness condition (2.15) for some e < ^. This 
constraint can, incidentally, be imposed on the system by choosing an appropriate 
gauge field action and one then obtains a lattice regularization of QCD which pre¬ 
serves the chiral and flavour symmetries without violating basic principles. As far as 
we know all other regularizations of QCD that have been considered to date break 
chiral symmetry and the folklore has been that this is in fact unavoidable. 

For numerical simulations of lattice QCD with the standard gauge action, the 
locality of Neuberger’s operator must be guaranteed for larger values of the gauge 
coupling as well. We have addressed this question in quenched QCD and did not 
find any indication that the locality is lost at the couplings of interest, provided the 
parameter s is chosen appropriately. 

In principle one may now use Neuberger’s operator to calculate the hadron spec¬ 
trum etc., but this may still be somewhat premature, because there are other choices 
of A which lead to significantly smaller localization ranges, at least in the free case 
[8,28]. This may be important in practice, since one cannot afford to simulate very 
large lattices. Moreover, as has been pointed out by Niedermayer [28], the proba¬ 
bility for near-zero modes may be very much suppressed for some of the proposed 
choices of A and this too could make the calculations easier. 

We are indebted to Peter Hasenfratz, Ferenc Niedermayer and Peter Weisz for 
encouragement and many helpful discussions. M. L. would like to thank the Institute 
for Theoretical Physics at the University of Bern for hospitality in June, where part 
of this paper has been completed. The computer resources for this project have been 
provided by CERN and CIEMAT. We thank these instutions for their support. 


Appendix A 


The notational conventions used in this paper are standard. We consider a four¬ 
dimensional hyper-cubic lattice with spacing a and variable size. If the lattice is 
finite we impose periodic boundary conditions although most results hold for other 
boundary conditions as well. The gauge field is represented by unitary matrices 
U {x, fi) where x runs through all lattice points and g = 0,... ,3 labels the space-time 
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directions. Dirac fields ^p{x) carry a Dirac and a colour index as in the continuum 
theory. The gauge covariant forward and backward difference operators act on such 
fields according to 


= ^[U (x, fi)^(x + afi) - V'(x)], (A.l) 

V*V'(x) = - [i){x) -U{x- afi, ny^'ipix - afi )], (A.2) 

a 

where jl denotes the unit vector in direction /j,. Since we are in euclidean space, the 
Dirac matrices can be taken to be hermitean, 

7m = (A-3) 

and our conventions for 75 and a^i, are 

75 =70717273, = 2 • (A.4) 

Repeated indices are always summed over unless stated otherwise. 


Appendix B 


The expansion in Legendre polynomials derived in sect. 2 is exponentially conver¬ 
gent. It is not obvious, however, that it may be differentiated with respect to the 
gauge field (or any other parameter), because the differentiated polynomials need not 
be uniformly bounded. In the following lines we establish a bound which excludes 
such an irregular behaviour. 

We first derive an integral representation for the Legendre polynomials. Starting 
from eqs. (2.2) and (2.3) we have 

Pk{z) = / ^ + 1), (B.l) 

where the integration runs along a circle in the complex plane centred at the origin. 
The radius r of the circle should be strictly less than 1 to avoid the singularities of 
the integrand. Because of the square root and since 2 ; is an operator, eq. (B.l) is 
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not easily differentiated. To overcome this difficulty we make use of a well-known 
identity to rewrite the integral in the form 

P,^{z)= [ — ([— -2wz + l + a^)-\ (B.2) 

J_oo tt J 2m 

Note that the denominator of the integrand can be factorized according to 

up' — 2wz -t- 1 + = (re — P){w — u), u = z + i[l + a^ — . (B.3) 

Since = 1 -|- it is then immediately clear that 

||(tc^ — 2 rt; 2 ;-|-1 -|-II < { (l-|-— r} (B.4) 

In particular, the integral (B.2) is norm convergent. 

Let us now assume that A'^A depends on some parameter r in a differentiable 
manner such that i = dzjdr has finite norm. An upper bound on the derivative 
of the Legendre polynomials with respect to r is then obtained by differentiating 
eq. (B.2) and applying eq. (B.4). The right-hand side of the resulting inequality 

\\Pkiz)\\<4\\z\\r~^ [ — {[l + - r] ^ (B.5) 

Jo ^ 

can be evaluated by substituting 

(1 -|- — 1 = (1 — r)p^, 0 < p < oo. (B.6) 

After some algebra one then ends up with the bound 

||-Pfc('2)|| < constant x ||i|| r“^(l — r)~‘^, (B.7) 

where the constant is independent of k and r. 

So far the radius r has not been specified apart from the requirement that it 
should be in the range 0 < r < 1. We may now adjust the radius so that the factor 
r“^(l — r)~^ is minimized. Using simple estimates this leads to the bound 

||7A(^)|| < constant x ||i|| (1 -|- k)'^. (B.8) 

The differentiated series (2.6) is hence exponentially convergent with the same ex¬ 
ponent as the original series. Similar estimations show that this is also true when 
higher-order differential operators are applied (each differentiation gives rise to an 
additional factor of (1 + kp in the bound on the Legendre polynomials). 
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Appendix C 


To establish eq. (2.16) we first expand the product A using simple identities for 
the covariant difference operators and the Dirac matrices. As a result one gets a 
sum of terms, 

AtA = l + + + (C.l) 

which can be treated separately (setting s = 0 has been essential here to ensure the 
cancellation of some non-trivial diagonal terms). Explicitly they are given by 

V^V^V,, (C.2) 

= ^ia^.a2[v; + V^,V*+V,], (C.3) 

= -7^a2 [V; + V^, V* - V,] . (C.4) 

The commutator terms (C.3) and (C.4) are proportional to the held strength and 


should thus be of order e. To prove this we note that 

[V^, V^] V'(x) = 

{C/(x, n)U{x + afi, v) — U{x, u)U{x + at>, ;u)}V’(a: + ajl + av). (C.5) 

The curly bracket in this equation is equal to a unitary matrix times 1 — U{p) for 
some plaquette p. Eq. (2.15) thus implies the bound 

||a'[V^,V,]|| <e (C.6) 

and the same inequality also holds if one or both forward difference operators are 
replaced by backward difference operators. In particular, 

IIC^i/ll < 2e and ||D^jy||<4e. (C.7) 

To bound the hrst term, eq. (C.2), we rewrite it in the form 

= a^V*- a^V* [V^, - V.], (C.8) 


which shows that it is equal to a non-negative operator plus another operator with 
norm less than 4e. Taken together these estimates imply that the right-hand side of 
eq. (C.l) is bounded from below by 1 — 30e which proves eq. (2.16). 
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